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Abstract 

The entanglement in quantum XY spin chains of arbitrary length is investigated via the geo- 
metric measure of entanglement. The emergence of entanglement is explained intuitively from the 
perspective of perturbations. The model is solved exactly and the energy spectrum is determined 
and analyzed in particular for the lowest two levels for both finite and infinite systems. The overlaps 
for these two levels are calculated analytically for arbitrary number of spins. The entanglement is 
hence obtained by maximizing over a single parameter. The corresponding ground-state entangle- 
ment surface is then determined over the entire phase diagram, and its behavior can be used to 
delineate the boundaries in the phase diagram. For example, the ficld-dcrivativc of the entangle- 
ment becomes singular along the critical line. The form of the divergence is derived analytically 
and it turns out to be dictated by the universality class controlling the quantum phase transition. 
The behavior of the entanglement near criticality can be understood via a scaling hypothesis, anal- 
ogous to that for free energies. The entanglement density vanishes along the so-called disorder 
line in the phase diagram, the ground space is doubly degenerate and spanned by two product 
states. The entanglement for the superposition of the lowest two states is also calculated. The 
exact value of the entanglement depends on the specific form of superposition. However, in the 
thermodynamic limit the entanglement density turns out to be independent of the superposition. 
This proves that the entanglement density is insensitive to whether the ground state is chosen 
to be the spontaneously Z2 symmetry broken one or not. The finite-size scaling of entanglement 
at critical points is also investigated from two different view points. First, the maximum in the 
field-derivative of the entanglement density is computed and fitted to a logarithmic dependence of 
the system size, thereby deducing the correlation length exponent for the Ising class using only the 
behavior of entanglement. Second, the entanglement density itself is shown to possess a correc- 
tion term inversely proportional to the system size, with the coefficient being universal (but with 
different values for the ground state and the first excited state, respectively). 

PACS numbers: 03.67.Mn, 03.65.Ud, 64.70.Tg, 05.70.Jk 



2 



I. INTRODUCTION 



Entanglement has been recognized in the past decade as a useful resource in quantum 
information processing jl-4|. Very recently, it has emerged as an actor on the nearby 
stage of quantum many-body physics, especially for systems that exhibit quantum phase 
transitions js-ll], where it can play the role of a diagnostic of quantum correlations. Quan- 



tum phase transitions 



12| are transitions between qualitatively distinct phases of quantum 



many-body systems, driven by quantum fluctuations. It is conventionally defined as the 
non-analyticity in the ground-state energy . In view of the connection between entangle- 
ment and quantum correlations, one anticipates that entanglement will furnish a dramatic 
signature of the quantum critical point as the ground-state wavefunction itself should also 
exhibit certain non-analyticity, which thereby manifests in the singularity of correlations. 
To describe the ground state usually requires exponentially many parameters, and therefore 
to detect the non-analyticity hidden in the wavefunction may require an ingenious design of 
observables. Ground-state energy is usually such a good indicator but it may fail to exhibit 



singularity, for example, near the isotropic antiferromagnetic point in the XXZ model |13| . 
where there is a continuous transition. Entanglement is thus expected to be an alternative, 
useful quantity in assisting the determination of quantum phase transitions, complementing 
traditional methods. On the other hand, the more entangled a state is, usually, though not 
necessarily, the more useful it is likely to be a resource for quantum information process- 
ing [l^. It is thus desirable to study and quantify the degree of entanglement near quantum 
phase transitions. By employing entanglement to diagnose many-body quantum states one 
may obtain fresh insight into the quantum many-body problem jsj. 

To date, progress in characterizing and quantifying entanglement has taken place greatly 
in the domain of bi-partite systems, and much understanding has be gained in the multi- 



partite settings, although a complete picture has not yet been achieved 



isj. Much of the 



previous work on entanglement in quantum phase transitions has been based on bi-partite 



measures, i.e., focus has been on entanglement either between pairs of parties 



m 



1611 or 



between a part and the remainder of a system [SMlOl. llTl . |18[ |. For multi-partite systems, 
however, the complete characterization of entanglement requires the consideration of multi- 
partite entanglement, for which a consensus measure has not yet emerged. 

Singular and scaling behavior of entanglement near quantum critical points was discov- 
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ered in important work by Osterloh and co-workers [7|, who invoked Wootters' bi-partite 
concurrence [lol as a measure of entanglement. The success of the use of concurrence has 
been understood via its connection to two-point correlation functions and the derivative of 



ground-state energy 



20l |. The drawback of concurrence is that it can deal with only two 



spins (each with spin-1/2) even though the system may contain an infinite number of spins. 
Although attempts have been made to generalize concurrence to many spin-1/2 systems via 
the time reversal operation, the generalized concurrence loses its connection to the entan- 
glement of formation 2l|. Furthermore, the use of concurrence may give rise to spurious 



transition points 



22|. 



Another approach is to consider the von Neumann entropy of a subsystem of L spins with 
the rest N — L spins of the system. It is found that for critical spin chains the entropy scales 
logarithmically with the subsystem size L for N ^ oo, with a prefactor that is related to the 
central charge of the corresponding conformal theory js iO]. The entanglement addressed 
in this case is only between a subsystem and the rest of the system, and the connection to 
central charge is interesting. 

A different approach without using entanglement is to consider fidelity or overlap between 
two different ground-state wavef unctions, which is now called fidelity measure. Non-analytic 
behavior in the fidelity measure is expected for system parameters close to critical points. 



23|. 



This approach has been shown to be successful 

Recently, Barnum and co-workers HI] have developed an entanglement measure, which 
they call generahzed entanglement. Instead of using subsystems they use different algebras 
and generalized coherent states to define the entanglement. They have also applied the 
generalized entanglement to systems exhibiting quantum phase transitions. Their approach 
opens a new approach to multi-partite entanglement. Nevertheless, there is no a priori 
choice of which algebra, amongst all possible ones, is the most natural one to use. 

In addition to spin chains other models that have been studied by using either the von 
Neumann entropy (also known as entanglement entropy) or the concurrence as the entangle- 
ment measure include: (i) the super-radiance model, in which many two-level atoms interact 



with a single-mode photon field 



24| : and (ii) the one-dimensional extended Hubbard model. 



in which electrons can hop between the nearest neighbors and there are Coulomb interactions 



among electrons on the same site and with nearest-neighbor electrons as well 



25| . Further- 



more, in addition to spin systems, entanglement entropy has been examined in both bosonic 
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and fermionic systems 18|, |26| and the area law (or the breaking of it) has been investigated. 



In view of other approaches that connect quantum information ideas to condensed matter 

n 

settings (or vice versa), Verstraete and co-workers [27] have defined an entanglement length, 
viz., the distance at which two sites can establish a pure-state entanglement at the cost 
of measuring all other sites. They found that this entanglement length is lower bounded 
the correlation length. All these, including the theme of the present paper, are aimed at 
approaching many-body problems from different, and hopefully fresh, perspectives, and at 
complementing traditional statistical-physical approaches. 

In the present paper, we expand and extend results of Ref. ^8\. We apply the global 



geometric measure that we have developed in Ref. [291], based on a geometric picture 
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33| : it provides a holistic, rather than bi-partite, characterization of the entanglement of 



quantum many-body systems. Our focus is on one-dimensional spin systems, specifically 
ones that are exactly solvable and exhibit quantum criticality. The geometric measure has 
been applied to many models in one dimension 34l-l38| as well as two dimensions 39|] and 
seems to become a useful tool as well. For these systems, it is found that entanglement 
behaves in a singular manner near the quantum critical points. This supports the view that 
entanglement — the non-factorization of wave functions — reflects quantum correlations. This 
will be explained in detail using the transverse-field XY models as an specific example. 

The remaining structure of the present paper is as follows. In Sec. HIl we introduce the 
geometric measure of entanglement (or in short, geometric entanglement) that we shall use 
to quantify the entanglement of a many-body state. In Sec. IIIIl we introduce the spin XY 
models in a transverse field, the focus of present paper. We explain in an elementary way the 
emergence of entanglement away from critical points, so as to expect the drastic behavior of 
entanglement near criticality. In Sec. llVt we diagonalize the Hamiltonian, invoking Jordan- 
Wigner and Bogoliubov transformations. We carefully explain the boundary conditions 
placed on the spin and hence the fermionic operators. In Sec. |Vl we derive the overlap 
of the lowest two levels with an Ansatz product state, the result of which is the basis for 
deriving entanglement for these two states, including the true ground state. Section I VI I 
contains the main results of the present paper. We present the results of ground-state 
entanglement over the entire phase diagram. In Sec. IVI At we illustrate the entanglement 
results for finite systems. In Sec. IVIBt we derive the expression of entanglement for infinite 
systems. There, the divergence of ground-state entanglement is analyzed and the connection 
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to the universality classes is established. In Sec. IVI C\ we show the vanishing of entanglement 
along the disorder line by demonstrating the ground space is doubly degenerate and spanned 
by two product states. In Sec. IVIDl we discuss the issue of entanglement superposition, and 
argue the irrelevance of spontaneous symmetry breaking in the determination of ground- 
state entanglement density. In Sec. IVI E\ we study the issue of finite-size effect to the 
entanglement density at criticality. Finally, we summarize in Sec. IVIIt and give a generic 
picture of how entanglement behaves near criticality, using a scaling hypothesis. How the 
geometric entanglement serves to detect non-analyticity of ground-state wavefunction is also 
discussed. 



II. GLOBAL GEOMETRIC MEASURE OF ENTANGLEMENT 



We quickly review the global measure that we shall use in the present paper. Consider 
a general, n-partite, normalized pure state: = X]pi---p„ '^piP2---pJ^pi^m ' ' '^v'n)- If the 
parties are all spin- 1/2 then each can be taken to have the basis {| t), | i)}- Our scheme 
for analyzing the entanglement involves considering how well an entangled state can be 
approximated by some unentangled (normalized) state (e.g., the state in which every spin 
points in a definite direction): |$) = ®'2=i l^*-*-*)- The proximity of to |$) is captured by 



their overlap; the entanglement of |\1') is revealed by the maximal overlap 29N33I. |40 | 



A^ax(^)=max|($|^)|; (1) 

the larger A^ax is, the less entangled is (Note that for a product state, Amax is unity.) If 
the entangled state consists of two separate entangled pairs of subsystems, Amax is the prod- 
uct of the maximal overlaps of the two. Hence, it makes sense to quantify the entanglement 



of \^) via the following extensive quantity 41| 



Eiog,(^) = -log2ALxW, (2) 

This normalizes to unity the entanglement of EPR-Bell and A/^-party GHZ states, as well 

as giving zero for unentangled states. Finite- iV entanglement is interesting in the context 
of quantum information processing. To characterize the properties of the quantum critical 
point we use the thermodynamic quantity £ defined by 

£ = lim £j^, (3a) 

N^oo 

£n = N-'Ei,^^{^), (3b) 
6 



where 8n is the entanglement density^ i.e., the entanglement per particle. 

For translation invariant state, it is shown that the entanglement is either completely 
in a product state or is globally entanglement 42|. Therefore, the above quantification of 
entanglement by comparing to product states makes a natural choice. 

We remark that by appropriate partitioning of |0) into various product forms, a hier- 



archy of entanglement can be obtained 
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recently also been applied to the XY model 



45 



44j . This hierarchical geometric measure has 
46| . One can generalize the form of product 



states, e.g., to a product of blocks of two spins (or multiple spins), and then calculate the 



entanglement per block 



36|. In fact, this approach enables the connection of the spatial 



renormalization group (RG) (i.e., coarse-graining) procedure to the geometric measure, and 



the entanglement under RG can be properly defined 



37|. 



III. QUANTUM XY SPIN CHAINS AND EMERGENCE OF ENTANGLEMENT 

We consider the family of models governed by the Hamiltonian 

^XY = - E (^-J-j^i + V-Ki + ^ ^l) > (4) 

where r measures the anisotropy between x and y couplings, h is the transverse external 
field, lying along the ^-direction, and we impose periodic boundary conditions, namely, a 
ring geometry. At r = we have the isotropic XY limit (also known as the XX model) 
and at r = 1, the Ising limit. All anisotropic XY models (0 < r < 1) belong to the same 
universality class, i.e., the Ising class, whereas the isotropic XX model belongs to a different 
universality class. XY models exhibit three phases (see Fig. H]): oscillatory, ferromagnetic 
and paramagnetic. In contrast to the paramagnetic phase, the first two are ordered phases, 
with the oscillatory phase being associated with a characteristic wavevector, refiecting the 
modulation of the spin correlation functions (see, e.g., Ref. 4^]). We shall see that the 
global entanglement detects the boundaries between these phases, and that the universality 
class dictates the behavior of entanglement near quantum phase transitions. 

Before we solve the entanglement of the XY model, we give perturbative analysis of, as 
an illustration of how entanglement arises and vanishes, the Ising model in a transverse field 



7 



[VIZ. r = i) 

N 



i=l 

At h = the ground state is that with all spins pointing up in the x-direction | — t-— >■ 
•■■ — or down | ■■■ which is manifestly unentangled. The ground state can 
be any superposition of (!—;■—;• ■■■ — and | ■■■ ^) when the Z2 symmetry is not 

spontaneously broken. For example, the states (| — ■■■ — )•) ± | ■■■ ^))/V2 are 
actually the two lowest levels obtained from solving the models using Jordan- Wigner and 
Bogoliubov transformations and they both have -Eiogj = 1- (For small h the entanglement 
rises quadratically in the case of unbroken symmetry instead of quartically, as we shall show 
shortly.) We shall see later that whether or not we use a broken-symmetry state has no 
effect in the thermodynamic limit. For small h (i.e., h ^ 1/\/N) one can obtain the ground 
state by treating the ha^ terms as perturbations. Take the ground state at /i = to be 
|— ;■— > ■ ■ ■ — )■). Then first-order perturbation theory for the ground state gives 

1 /, . h 



1 I Nh'^ 
16 



■ + (6) 



Using the method described in Ref. 



29| we obtain -Eiogj ^ Nh"^ / {16\ia2) to leading order 



in h. At /i = 00 the ground state is a quantum paramagnet with all spins aligning along 
the external field: Iff ... f), and once more is unentangled. To 0{l/h) perturbation theory 
gives (treating erf a^'^.i terms small) 

=^(iit-t) + ^Eit-w«...t>). (7) 

for which -E'logj ^ A^/(16/i^ ln2), to leading order of 1/h. The quantum phase transition 
from a ferromagnetic to a paramagnetic phase occurs ai h = 1 12|. The two lowest levels, 
which we denote by |\I'i/2) and |\E'o) (for reasons to be explained later) are, respectively, the 
ground and first excited states, and they are asymptotically degenerate for < h < 1 when 
N ^ 00. 



IV. DIAGONALIZATION OF THE HAMILTONIAN 



As is well known 12|, |47|, |49|, the energy eigenproblem for the XY spin chain can be 
solved via a Jordan- Wigner transformation, through which the spin degrees of freedom are 
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recast as fermionic ones, followed by a Bogoliubov transformation, which diagonalizes the 
resulting quadratic Hamiltonian. 

The Jordan- Wigner transformation that we shall make from spins (ex's) to fermion par- 
ticles (c's) is 

af = l-24Q, (8a) 

i-l 

^! = -^Ili^-^4cj){c,-cl). (8c) 

One has to pay attention to the boundary conditions that are to be imposed on the c's. 
Although periodic in the cr's, one cannot simply take 

AT 

<+l = n (1 ~ Cj) (Cjv+l + = (7i = (ci + c{) , (9) 

and conclude that either (i) HjLi (l~2c] Cj) = 1 and cjv+i = Ci, or (ii) Y[f=i (^^^c] Cj) = — 1 
and cat+i = — Ci, neither of which are correct if one wishes to obtain the correct spectrum 
and eigenstates for arbitrary finite N. Instead one should impose that (as the (A^ + l)-th 
site is identified as the first site) 

which then leads to 

N 

{CN + c!v) (cAT+l + c!v+i) = - n (1 - 2c] Cj) {cn + cjv) (ci + c\) . (11) 

The two possible conditions that satisfy this are either (I) Y[f=i (l ~ ^c] cj) = — 1 and 
cn+1 — ci, or (II) Y[f=i ~ 2c] Cj) = 1 and cn+i = — ci. The operator 

N 

]^(l-2c]c,) =e'"^^-"^^- (12) 

counts whether the total number of particles is even (+1) or odd (—1). For c's that are 
periodic, the number is odd, whereas for antiperiodic c's, this number is even. 
To incorporate these two boundary conditions on the c's, we take 

N-l 



Cj 

m=0 



1 (13) 



0.5 1 1.5 2 2.5 

h 



FIG. 1: Energy spectrum vs. magnetic field h for the Ising spin chain with N = 4 qubits from the 
Bogohubov diagonahzation. The results are identical to exact matrix diagonalization. 

where 6 = for periodic c's; 6 = 1/2 for anti-periodic c's. (This explains why we label the 
lowest two states by 1^1/2) and |^o)-) The momentum index m ranges from to A?" — 1. In 
terms of these fermion operators the Hamiltonian becomes 

N-l 



n = -Nh-J2 \ [2cos-^(m + 6) -2/i]( 

m=0 

+ir sin f (m + 6) [cg^c^^l^_,, + } ■ (14) 

Upon using the Bogoliubov transformation 

c(^) = cos e^^V^^ + t sin ^^)7S^i„.-2. , (15) 

with 

ta.29S'=r.in?l%±^/(fe-co.B^(!^), (16) 



one arrives at the diagonal the Hamiltonian: 

N-l 



n^-Nh + J2et^{^(^^(^^-l), (17a) 



m=0 



g^^ = 2y(/^-cos +r^sm'^-—^, (17b) 

except for the special case that £q°^ = 2(h — 1). 

We remark that we have not left out any constant in diagonalizing the Hamiltonian in 
either case, so the energy spectrum is exact. For each value of b the diagonahzation gives 2^ 
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0.5 1 1.5 2 2.5 

h 



FIG. 2: (Color online) Energy gap vs. magnetic field h for the Ising spin chain with A'^ = 4 (red 
shorter dashed), = 10 (blue longer dashed), and = 100 (black solid curve). 

energy eigenvalues, so there are 2^"*"^ in total. Half of them are spurious. In determining the 
correct 2^ states from the 2^+^ solutions, one has to impose a constraint from the boundary 
conditions. Namely, in case I there can be only odd number of fermions, whereas in case II 
there can be only even number of fermions. 

For 6 = 0, viz, the odd-number-fermion case, the lowest state |\E'o) is such that 
(7m^^7m^) = except that (7o°''^7o°'') = 1- Its energy eigenvalue is 

Ei;'\r,h) = {h-1)- J2\l(h- cos ^) +r^sm'^. (18) 

m=l 

For b = 1/2, namely, the even-number-fermion case, the lowest state |^I^i/2) is such that 
(7m^^^^7m^^^) = for all m. Its eigen-energy is 

m=0 

We see that, as — t- oo, the above two energy levels are degenerate for h <1. Further- 
more, as A^ — )■ oo the difference between the two energy levels becomes 

Ef\r, h) - E^^/'\t, h) = 2{h - l)Q{h - 1), (20) 

where 6(x) = 1 if x > and zero otherwise. The way the energy gap vanishes as /i — )■ 1"*" 
gives a relation between two exponents 

zu = 1; (21) 
11 
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0.5 1 1.5 2 

h 

FIG. 3: Energy Energy difference dE = E^^^ — E^^'^^ vs. magnetic field h for the XY spin chain 
with r = 0.2 and = 8 qubits from the Bogohubov diagonahzation. 

z is the dynamical exponent (defined via the vanishing of energy gap A ~ \h — h^Y''') and u 
is the correlation- length exponent (defined via ~ \h — /id"^)- 

We show in Fig. [T]the energy spectrum for = 4 Ising spin chain obtained from the above 
JW and Bogoliubov diagonahzation, which is identical to direct matrix diagonahzation (not 
shown). This means that the choice of our boundary conditions below Eq. ffTTj) is correct. 
We illustrate the gap in the spectrum for the Ising chain in Fig. [21 which is obviously 
degenerate for h < 1 and linear {h — 1) for h > 1, as derived in Eq. ([20]). We see that E^^ 
is asymptotically (i.e., as — ?■ oo) degenerate with Eq^^"^^ in h E [0, 1], but the former has 
higher energy than the latter for h > 1. For finite A^, i?Q°^ is larger than E^^'^\ However, 
this is not the case for general XY model, as we illustrate in Fig. [3] for r = 0.2, where we see 
that the two branches switch roles of being the true ground state, depending on the value 
oih. 

Even though the ground state may switch back and forth between |^I/i/2) and |\l/o) for 
/i < 1, for /i > 1 |\I'i/2) is the ground state. In particular, at /i = 1 we calculate AAr(r) = 
E^\r^ h = 1) — E'^^'^\r^ h = 1) for the finite chain with A^ spins and find that 

A^{r) = ^ + 0{1/N'). (22) 

It is interesting to see that the difference in the two energies is proportional to the anisotropy 
r and, as expected, inversely proportional to the system size A^. 
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V. DERIVATION OF OVERLAP OF THE LOWEST TWO STATES WITH THE 
PRODUCT ANSATZ STATE 



Having found the lowest two eigenstates, the quantity A^ax of Eq. ([T]) — and hence the 
entanglement — can be found, at least in principle. To do this, we parametrize the separable 
states via 

N r -J. 

1$) = ® [coste/2)|t)i + e'^»sin(e./2)|;),], (23) 
1=1 

where | t/i) denote spin states parallel/antiparallel to the 2;-axis. Instead of maximizing 
the overlap with respect to the 2N real parameters {^j, for the lowest two states it is 
adequate to appeal to the translational symmetry of and the reality of the ground-state 
wavefunctions. Thus taking C,i = ^ and (pi = we make the Ansatz: 

|$(0) = e-'i ^^-1 m • • • t) (24) 



for searching for the maximal the overlap Amax(^) (48|]. This form shows that this separable 
state can be constructed as a global rotation of the ground state at /i = oo, viz., the separable 
state I 'l^'l" . . . 'I") . In this particular limit the boundary condition on the c's is irrelevant, as 
the dominant term in the Hamiltonian is — h{l — 2cjcj). 

The energy eigenstates are readily expressed in terms of the Jordan- Wigner fermion 
operators, and so too is the family of the Ansatz states |<I>(^)). By working in this fermion 
basis we are able to evaluate the overlaps between the two lowest states and the Ansatz 
states, presented in the next section. 

We first analyze 6=1/2 case, viz., the even-fermion case. The lowest state \^ii2{r,h)) 
has zero number of Bogoliubov fermions. It is related to the state that has no c-fermions, 
i.e., = 1 1 ■ ■ • t) via 

q(l/2). ,^ -(l/2)t-(l/2)t 



\^,/,{r,h)) = H cos^^/2)(r,/i)e^*^^^-"(^''^)^- "'^--^-11]) (25a) 

m=0 

Y[ ^cos 9m{r, h) +i sin 9m{r,h)cl^c^^_^_-^ (25b) 



m=0 
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The Ansatz state is then 



\m) = n [^^4 + '^"^ I n (1 - 2c;q)(ci - c,)] i^) (26a) 

j=i i<i<j 

= IK^o^+'^'^r^^^'^^ ^^^^^ 

= cos^ ^ e*^° 24... (26c) 

where we have suppressed the index (1/2). The term 'YlijKi can be rewritten in momen- 
tum space as 

E M = i E -t^^^^44_. (27) 

l<j<l<N m=0 

Thus, for even N 

|$(0) = (l + tan|5^ct) n cos^l + zsin^cot^^i^^cLgjv^Jl^^ (28a) 

j=l m=0 ^ ^ 

whereas for odd 



|$(e))=(l + tan|^cj.)cos| J] cos^ ^ + z sin^ | cot ^^^^44_i ) ^ 

j=l m=0 ^ ^ 

(28b) 

Therefore, the overlap of the state \'ifi/2{r, h)) with |$(C)) for even N is 



{^y,{TMm)= n fcos^(i/2)(r,/i)cos2^ + sin^^/2)(r,/i)sin2^cot^^^^yJ^ 

m=0 ^ 

(29a) 

whereas for odd N 

(*V2(r, h)\m) = COS I n (co^^^^'Hr, h) cos' | + sinC/2)(r, h) sin^ | cot !Li!!!±l2 

m=0 ^ 

(29b) 

Next, we discuss the 6 = (odd-fermion) case. The lowest allowed state is the one with 
one 7q°^ = Cq°^ fermion: 

\Mr,h)) ^ 7riG(r,/i)) = c^^^^\Gir,h)), (30) 
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where \G{r, h)) is the state with no 7 fermions: 

m<f 

\G{rM - n [cos^W(r,/i) + isin^W(r,/i)cWtg(J)_g|0). (31) 



m=l 

Similar to the h — 1/2 case, by using 



m<f 



E c]c\ = ^Y.^ot^-^^''cfl. (32) 

l<j<l<N m=l 



we obtain that for even A?" 



I*(0)=(l + ^tan|4)cos^^ n +^sin'fcot^44-^)l^^), (33a) 

m=l ^ ^ 



whereas for odd N 



l*(e))= (l + V^tan|ct)cos| J] ( ^o^' f + ^ ^ ^44-^ ) 1^^)- (33b) 

m=l ^ ^ 



Therefore, the overlap of \^o{i^: h)) with ^or even is 

(*o(r,/i)|$(0) = v^sin| n (cose(r,/i)cos^| + sin^W(r,/^)sin^|cot^], (34) 

m=l ^ ^ 

whereas for odd N 

(*o(r,/i)|$(0) = v^sin^cos^ J] (cos^i^H^,^) cos^ | + sin ^ (r, /.) sin^f cot^J . 

m=l ^ ^ 

(35) 

Summarizing the above derivations: with |^o) (|^i/2)) denoting the lowest state in the 
odd (even) fermion-number sector, we arrive at the overlaps 

{Mr,h)\m)-fS\0 n [cose^)(r,/i)cos^(e/2) + sin^W(^,/,)si„2(^/2)cot(A;Sy2)] , 

m=l-2& 

(36) 



with 



^ '^-^(rn + b), tan2^£)(r,/.) ^ rsmk^^y{h-cosk^j^]^y, (37a) 
fi^/'\0 = 1, fS\0 = VNsm{^/2) cos(e/2), (N even); (37b) 
/S/^'(0=cos(e/2), fS\0 = y^MC/2), (AT odd); (37c) 
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where b = 0, 1/2 and m G [0, — 1] is the (integer) momentum index. The above results 
are exact for arbitrary A^, obtained with periodic boundary conditions on spins rather than 
in the so-called c-cyclic approximation 49|. Given these overlaps, we can readily obtain 
the entanglement of the ground state, the first excited state, and any linear superposition, 
cosa|\E'o) + sina|\E'i) of the two lowest states, for arbitrary (r, h) and A^, by maximizing the 
magnitude of the overlap with respect to the single, real parameter ^. 

We remark that with the lowest states explicitly expressed in terms of the fermionic 
language, their fidelity measure |23| can be easily evaluated. Our focus here is the geometric 
measure of entanglement, which we derive in the next sections. 



VI. ENTANGLEMENT IN THE TRANSVERSE-FIELD XY CHAINS 

The formulas regarding the overlaps of |^E'6)'s with product state [iii Eqs. (l36l) 

and ( 1371) ] contain all the results that we shall explore shortly, including both finite-size 
and infinite-size limits. By analyzing the structure of Eq. (15B|) . we find that the global 
entanglement does provide information on the phase structure and critical properties of 
the quantum spin chains. Two of features, as captured in Figs. HI and 15], are: (i) although 
the entanglement itself is, generically, not maximized at the quantum critical line in the 
(r, h) plane, the field- derivative of the entanglement diverges as the critical line h = 1 is 
approached; and (ii) the entanglement vanishes at the disorder line r"^ + h"^ = 1, which 
separates the oscillatory and ferromagnetic phases. 

As is to be expected, at finite the two lowest states |\E'o) (with energy -Eq^^) and |\E'i/2) 
(with energy e'^^'^^) featuring in Eq. (l36l) do not spontaneously break the Z2 symmetry. 
However, in the thermodynamic limit they are degenerate for /?. < 1, and linear combinations 
are also ground states. The question then arises as to whether linear combinations that 
explicitly break Z2 symmetry, i.e., the physically relevant states with finite spontaneous 
magnetization, show the same entanglement properties. In fact, we see from Eq. fl36l) that, 
in the thermodynamic limit, overlaps for |\l/o) and |^E'i/2) are identical, up to the prefactors 
and f^^'^^. These prefactors do not contribute to the entanglement density, and the 
entanglement density is therefore the same for both |\E'o) and |^E'i/2)- It further follows that, 
in the thermodynamic limit, the results for the entanglement density are insensitive to the 
replacement of a symmetric ground state by a broken-symmetry one. 
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FIG. 4: (Color online) Entanglement density (upper) and phase diagram (lower) vs. {r,h) for the 
XY model with N = 10^ spins, which is essentially in the thermodynamic limit. There are three 
phases: O: ordered oscillatory, for + /i^ < 1 and r 7^ 0; F: ordered ferromagnetic, between 
+ /i^ > 1 and /i < 1; P: paramagnetic, for h > 1. As is apparent, there is a sharp rise in the 
entanglement across the line h = 1, which signifies a quantum phase transition. The arc 
along which the entanglement density is zero (see also Fig. [5]), separates phases O and F. Along 
r = lies the XX model, which belongs to a different universality class from the anisotropic XY 
model. 17 
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FIG. 5: (Color online) Entanglement density and its /i-derivative (inset) for the ground state of 
three systems at = oo. Solid line: Ising (r = 1) limit; dashed line: anisotropic (r = 1/2) XY 
model; dash-dotted line: (r = 0) XX model. For the sake of clarity, the XY-case curves are shifted 
to the right by 0.5, indicated by the arrow. For the r = 1/2 case, at the entanglement density 
vanishes at h = \/l — r'^, which is a general property for the anisotropic XY model. Note that 
whilst the entanglement itself has a nonsingular maximum at h ~ 1.1 (Ising), h ~ 1.04 (r = 1/2 
XY), and h = (XX), respectively, it has a singularity at the quantum critical point at /i = 1, as 
revealed by the divergence of its derivative. 

One the other hand, entanglement of finite-size systems does depend on the exact ground 
state(s) and which superposition to take when there is degeneracy. For instance, the ground 
state switches between |\[^o) and |\I/i/2) (see Fig. [3]), and the ground-state entanglement for 
a finite system can be broken into pieces (see Fig. [8]). Furthermore, we shall also investigate 
the entanglement of the superposition: cos^|\E'i/2) + sin6'|\l/o). 

A. Entanglement for finite spins 

Before we discuss the thermodynamic limit of the entanglement density, we compare 
the entanglement obtained via the results in Eq. fl36l) and that obtained via numerically 
diagonalizing the Hamiltonian and calculating the maximal overlap. In Figure [6] the results 
via each method are shown for the Ising case (r = 1) with small numbers of spins (A^ = 13 
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FIG. 6: Total entanglement vs. magnetic field h for ground (upper panel) and first excited (lower 
panel) states with small N = 13, 16, 19, 22 (from bottom to top) for the transverse Ising model 
(r = 1). The numerical results are shown as discrete points whereas the analytical results are 
shown as lines. This demonstrates that the analytical results are exact, even for small A^, both 
even and odd. 

through 22), it is seen that our analytical results are exact even for small A^, both even and 
odd. 

As another example, we show in Fig. [7| the total entanglement for the ground state of the 
XX chain with N = 100 spins. The entanglement varies with h stepwise. This makes sense 
as at /i = 0, the ground state has continuous symmetry in XY plane, and hence has high 
entanglement. As h increases, the z-component total angular momentum increases by one 
once the increase exceeds certain threshold, breaking the continuous symmetry, and hence 
the entanglement decreases stepwise until at /i = 1 when all the spins are pointing in the 
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FIG. 7: Total entanglement vs. magnetic field h for the ground state of XX chain with A'^ = 100 
qubits. 

z-direction. In the limit N ^ oo the change in the step per spin becomes infinitesimal and 
the curve becomes smooth shown in Fig. [S] as expected. 

As a final example in this section, we illustrate the strange feature of piecewise continuity 
for finite-system ground-state entanglement. Take r = 0.2 and = 8 for example. The 
ground state switch between the two states labeled by the energy Eq^^ (odd fermion number) 
and Eq^^"^^ (even fermion number), as already shown in Fig. [31 In Fig. [HI we show the 
entanglement for these two states (solid red and dashed blue dashed curves), as well as the 
ground-state entanglement obtained from exact diagonalization (shown as dots). The ground 
state for this system has a piecewise continuity in its entanglement. A similar behavior 
has also been observed in the context of the hierarchical geometric entanglement 45|, [46 1. 



We remark that when one considers the entanglement per site, the difference between the 
entanglement density of the two states vanishes in the thermodynamic limit, as we shall see 
below. 

B. Entanglement density of infinite chains 

From Eq. ( l36l) it follows that the thermodynamic limit of the entanglement density is 
given by 

2 

£{r, h) = —- — max / In [cos 9{fi, r, h) cos^(^/2) + sin 6'(yU, r, h) sin^(^/2) cot vr/i] ,(38) 
In 2 c Jo 
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FIG. 8: (Color online) Total entanglement vs. magnetic field h of the XY chain with r = 0.2 and 
N = 8 for the ground state (black dots), the minimum energy state in odd-fermion sector E^^^^ 
(dashed blue curve) and even-fermion sector E^^'^'' (sold red). We see that the ground state, and 
hence its entanglement, switches between these two states; see also Fig. [3l 

where tan2^(yU, r,h) = r sin 27r/i/ {h — cos 2n fi). 

Figure [5] shows the thermodynamic limit of the entanglement density S{r,h) and its h- 
derivative in the ground state, as a function of h for three values of r, i.e., three slices 
through the surface shown in Fig. HI As the r = 1 slice shows, in the Ising limit the 
entanglement density is small for both small and large h. It increases with h from zero, 
monotonically, albeit very slowly for small h, then swiftly rising to a maximum at h ^ 
1.13 before decreasing monotonically upon further increase of h, asymptotically to zero. 
The entanglement maximum does not occur at the quantum critical point. However, the 
derivative of the entanglement with respect to h does diverge at the critical point h = 1, as 
shown in the inset. The slice at r = 1/2 (for clarity, shifted half a unit to the right) shows 
qualitatively similar behavior, except that it is finite (although small) at h = 0, and starts 
out by decreasing to a shallow minimum of zero at h = \/\ — r'^. By contrast, the slice at 
r = (XX) starts out at = at a maximum value of 1 — 27(7/(71 In 2) ^ 0.159. (where 
7(7 ~ 0.9160 is the Catalan constant), the globally maximal value of the entanglement over 
the entire (r, h) plane. For larger h it falls monotonically until it vanishes aih = 1, remaining 
zero for larger h. 

Let us analyze the behavior of the entanglement derive for r 7^ (Ising universality class) 
case first. 
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1. Divergence of entanglement-derivative for the anisotropic XY models 



The starting point is Eq. ( 138|) . in which there is a maximization over the variable ^. The 
function to be maximized is 

F{^,r,h) = / ' rf/i In [cos ^(/i, r, /i) cos^(^/2) + sin ^(/i,r, /i) sin2(^/2) cot vr/i] . (39) 
Jo 

To find the stationarity condition, we demand the derivative with respect to ^ vanishes 
0), where 



{d^F{^,r,h) 



1 . ^ n d cos6{fi,r, h) — sm6{fi,r, h) cotiTfi 

5=5* 2 Jq cos6'(/i,r, /i) cos^(^/2) + sin6'(/i,r, /i) sin^(^/2) cotvr/i 



5=C* 
(40) 



Denote by ^*(/i) the solution for fixed r. Then the field-derivative of the entanglement is 

-2 



d,£ir,h) = ^d,FiCih),h) = ^ 



^^d^F{^,h)\^^+d,FiC,h) 



In 2 



dhFiCh), 



(41) 

where the first term in the square bracket vanishes identically due to the condition (l40l) . 
Thus (dropping the * on ,^ for convenience), 
2 



dh£ (r, h) 



In 2 

2 

'hr2 



dhFiC,h) (42) 
dh cos 9{fi, r, h) cos^(^/2) + sin 6{fi, r, /i) sin^(,^/2) cot iTfi 



cos 6'(/i, r, h) cos2(^/2) + sin 6'(/x, r, /i) sin^(^/2) cot vr/x 



Recall that tan 26'(/i, r,h) = r sin 27r/i/ (/i — cos 27r/i) and thus 



cos = v^(l + cos2^)/2, sin ^ = ^/{l -cos2e)/2, 

/i — cos 27ryU 



cos 26'(/i, r, /;,) 



^y (r sin 2tthY + {h — cos 2Tc^y 



-■ (43) 

(44a) 
(44b) 



Putting everything in Eq. ( 143|) . we get 



4^ (r, /;,) 



In 2 



sin 27r/i 



(rsin27r/i)^ + (/i — cos27r/i)^ 



V7 



(/i - cos 27r/i) cos^ (^ /2) - W V 



(/i — cos 27r/i) sin^ (^ / 2) cot vr 



45) 



^J \/ '+(/;, — cos 27r/i) cos^(,^/2) + y^V " — (/i — cos27r/i) sin^(,^/2) cot tt/x 

where \/ = a/ (?" sin + {h — cos 2'7r/i)2. 

We aim to explore the behavior near h = 1. First consider /i > 1 and define e = h — 1, 
which is the deviation from the critical point. Make the change of variables t = h — cos27rz/. 
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giving lower and upper limits e and 2 + e, respectively. We further shift the integration 
variable by e, arriving at 

1 



dhE{r,h) - 27rln2/ - r^)^^ + 2{r^ + e)t + 

Vi^ V -{t + e) cos2(^/2) - +{tVe)sm\^/2)V2^ 

-, (46) 



Vi^V +(i + e) cos2(^/2) + - (t+^sin2(e/2)V2^' 

where V = ^(1 - r2)t2 + 2(r2 + e)i + e^. 

As we inspect the limit /i — > 1 or e — > 0, we see that the above expression diverges, 
with the contribution coming from t small, i.e., infrared divergence. Large t{< 2) does 
not contribute to the divergence. Note further that only the second term in the numerator 

contributes to the divergence. We then proceed to evaluate the integral by separating it into 
two parts: 

/»2 /»5 /»2 

JO ^0 Js 

with 5 <^ 1. In the first region we only need to keep t to first order at most. Also noting 
that for t,e <^ 1, the first term in the denominator is much smaller than the second term, 
we get 



27rln2 



dt , ^ / =-. (48) 

2r2(t + j^) ^^2{r^ + e)t + -{t + e) 



Next, we simplify the second term (ignoring e when there is no danger in doing so) 



2(r2 + e)t + e2 + (i + e) Jt+^ t + e 



^J^2{r^ + e)tT^-{t + e) ^ 

Observing that only the first term on the right-hand side contributes to the divergence, we 
have that the divergent part is 



r -1 v^ + i^ 1 r^'"''' 1 1 

dt , V J"" ^ -— dt^^^. (50) 

27rln2 7o 2r2(i + ^) ^ft 4r7rln2 7o ^/t+l^ft ^ ' 

The divergence part is then (for 52r^/e^ ^ 1) 

1 1 , 1 / , \ 1 In (2^52^) , ^ 

- — — 2sinh-V^2r2/e2 ^ - — —In [2J52ryeA = -- — In e + — ^— — 51 

4r7rln2 ^ ' 2r7rln2 V ^ ' / 2r7rln2 2r7rln2 ^ ' 
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As the integral (H6|) does not depend on the choice of S, the part that involves S must be 
cancelled by the second half of the integration J^, which can be verified by direct evaluation. 
Therefore, for h very close to he = 1, "we have 

d£ 1 

-^\n\h-l\. (52) 

dh 27rrln2 ' ' ^ ^ 

In deriving the above divergence form, we have assumed that r ^ 0. Similar consideration 
can be applied to the case when h approaches 1 from below, and the behavior is the same. 



2. Divergence of entanglement-derivative for the XX limit of the model 



We now analyze the r = isotropic (XX universality class) case. It turns out that the 
analysis for this case is much simpler. To see this, we first note when r = we have the 
simplification 

cos29{fi, h) = sgn(/i — cos27r/i). (53) 

The above expression changes sign when h = cos27rz/. So let us introduce the variable 
fj,Q = (cos^^ h)/{2-n). The expression for the entanglement density Eq. (l38i) becomes 

2 

- max 
' ' '0 



In 2 ? 
2 



S^)—- — max / d/i In [cos 6'(/i, r, /i) cos^(^/2) + sin 6'(/i, r, /i) sin^(,^/2) cot TT/i] (54a) 



max 



In 2 5 



MO 



(i/i In sin^(,^/2) cot tt/x 



dfi lncos^(^/2) 



'0 ^MO 

Demanding stationarity with respect to C, gives the condition 



[f^o - ^(1 - cos ^)] sin ^ = 0. 



(54b) 



(55) 



The solution ^ = gives the entanglement density for h > 1. It is straightforward to see from 
Eq. fl54al) that the entanglement density is identically zero. The solution /iQ — 1(1 — cos^) = 
gives 



cos^ = 1 cos ^ h. 



TT 



This in turn gives the entanglement density for < /i < 1: 



S{h) 



In 2 



+ - ln(l — 2/xo(/i)) + / d/i In cot vr/x 
2 Jo 



(56) 



(57) 



2/io(/i) ^ JO 

Thus we have derived the entanglement density as a function of the magnetic field h in the 
XX limit. The result is shown in Fig. ([5]). 
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FIG. 9: Finite-size scaling. d8N/dh\h^^^ ^ vs. IniV, for N = 10"^, 3 x 10'', 5 x 10^,8 x 10'^, and 10^ 
(points) with r = 0.1. The sohd Hne represents the fit OSn /dh\h^^ ~ 2.30 In — 6.95. 

We see from the figure that the entanglement density at /i = is the highest, which 
being £{0) = 1 — 27c/(7rln2) ^ 0.159 by evaluating Eq. flFTI) at h = 0. The constant 
7c ~ 0.9160 is the Catalan constant. The entanglement density decreases monotonically as 
h increases until h = 1 beyond which it becomes zero identically. This qualitative behavior 
can be understood as follows. As the total 2;-component spin is conserved, increasing h 
simply increases the z-component spin of the ground state until h = 1 where all the spins 
are aligned with the field, hence there is no entanglement beyond this value of h. 

By directly taking the derivative with respect to h, we get 



dhS{h) = ^ In 

^ ^ 7rln2v/r^ 



IT — cos~^ h\ 1 — h 



cos ^ h 1 + h 



(58) 



Near h ^ 1, we have (putting 1 + h = 2 and evaluating the limit in the argument of log 
function) 

This completes our derivation of the singular behavior of the entanglement density near the 
critical points. 



3. Singularity of entanglement and quantum criticality 

Let us summarize the above derivations of the singular behavior in the entanglement 
density and discuss the connection to quantum criticality, in particular the correlation length 
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critical exponent. 

(1) The singular behavior of the field-derivative of the entanglement density ( 138|) in the 
vicinity of the quantum critical line acts as (for r 7^ 0) 

BE 1 

From the arbitrary-A^ results (136|) of the entanglement we analyze the approach to the 
thermodynamic limit, in order to develop further connections with quantum criticality. We 
focus on the exponent z/, which governs the divergence at criticality of the correlation length: 
Lc ^ \h — . To do this, we compare the divergence of the slope dE^ jdh (i) near h = 1 
(at = 00), given above, and (ii) for large at the value of h for which the slope is 
maximal (viz. /imax,iv), i-e., OEn /dh\h^^^ ^ 0.230r~MnA^ + const., obtained by analyzing 
Eq. fl5Bl) for various values of r; see Fig. [Hlfor the example of r = 0.1 case. Then, noting that 
(27rln2)~^ ^ 0.2296 and that the logarithmic scaling hypothesis 5^ (or see Appendix Ell 
identifies p with the ratio of the amphtudes of these divergences, 0.2296/0.230 ~ 1, we 
recover the known result that u = 1. Moreover, from Eq. (|2T1) we can extract the value of 
the dynamical exponent: z = 1. 

(2) Compared with r 7^ case, the nature of the divergence of dE/dh at r = belongs 
to a different (XX) universality class: 

From this divergence, the scaling hypothesis, and the assumption that the entanglement 
density is intensive, we can infer the known result [1^ that the critical exponent z/ = 1/2 for 
the XX model. Moreover, from Eq. (1211) we can extract the value of the dynamical exponent 
z = 2 for the XX model. 

In keeping with the critical features of the XY-model phase diagram, for any small but 
nonzero value of the anisotropy, the critical divergence of the entanglement derivative is 
governed by Ising-type behavior. It is only at the r = point that the critical behavior 
of the entanglement is governed by the XX universality class. For small r, XX behavior 
ultimately crosses over to Ising behavior. 
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C. Vanishing of entanglement density along the disorder line 



We find that along the hne r"^ + h"^ = 1 the entanglement density vanishes in the ther- 
modynamic limit. In fact, this line exactly corresponds to the boundary separating the 
oscillatory and ferromagnetic phases; the boundary can be characterized by a set of ground 
states with total entanglement of order unity, and thus of zero entanglement density. The 
entanglement density is also able to track the phase boundary {h = 1) between the or- 
dered and disordered phases. Associated with the quantum fluctuations accompanying the 
transition, the entanglement density shows a drastic variation across the boundary and the 
field-derivative diverges all along h = 1. The two boundaries separating the three phases 
coalesce at (r, h) = (0, 1), i.e., the XX critical point. Figures H] and reveal all these features. 

One extreme limit is the Ising case, i.e., r = 1 and h = 0, where the ground state is either 
|— )■ ■ ■ ■ — )■) or |-(— ^ ■ ■ • ^), both of these being unentangled. Any superposition of them 
is also a valid ground state, but it has entanglement of order unity. In the thermodynamic 
limit the entanglement per spin is identically zero. Is this a general feature along the disorder 
line? Before we establish this, recall that the energies of the lowest two levels are given in 
Eqs. (1181) and (1191) . Evaluating them at = 1 — /i^, we immediately find that both are — A^. 

Now let us evaluate the expectation value of the Hamiltonian (jl]) with respect to a 
separable state with all spins pointing in the same direction: 

{H) = -N (^(^^)^ + ^-^W'f + hW')) ■ (62) 

Denoting x = (cr^), y = (o"^), z = {cr^), we find that the above expression achieves its 
minimum value —N at r"^ + h"^ = 1 when 



x,y,z)= {±J-—,0,J-—]. (63) 




Therefore, the separable state satisfying the above conditions is the ground state, and there is 
a degeneracy. Hence, along the disorder line the total entanglement can be at most of order 
unity (due to superposition of the above degenerate states) and hence the entanglement 
density vanishes. We note that the theory of factorized ground states has recently been 



extensively developed and greatly understood [l6|, ISjJ, |52| , and it turns out that they occur 



quite often and their occurrence is related to the existence of quantum critical points nearby 
in the system parameter. 
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FIG. 10: (Color online) Total entanglement (E'logj) for an Ising chain vs. the field h for A'^ = 10 
and various superposition of cos0|^'i/2) +sin0|^'o): ^ = (solid black), 6 = vr/S (short-dashed 
red), 6* = 7r/4 (long-dashed blue), and 9 = it/2 (mixed-dashed green). 

D. Entanglement of superposition 

In this section, we discuss the entanglement of superposition for the two lowest states 
1^1/2) and 1^0): 

|^(r,/i,e)) = cose\'^i/2ir,h)) + sme\^oir,h)). (64) 

Why does this interest us? Take for instance the Ising model, r = 1. The states |^E'i/2) and 
1^0)) respectively, analytically continues to (as h 0) 

1*1/2(1, 0)) = (I >) + \^< ^))/v^ (65a) 

1*0(1, 0)) = (I . . . ^) - I . . . ^))/V2. (65b) 

This is a result of unbroken Z2 symmetry and both states possess global geometric entan- 
glement of unity. In the large system size limit, this symmetry is usually broken due to 
fluctuations of a random magnetic field, and the resulting symmetry-broken states will be 
unentangled, having zero entanglement. These symmetry-broken states can be represented 
as certain superposition of the original symmetry-unbroken states. But to evaluate the en- 
tanglement for these symmetry-broken states at a different field value h, we need to be able 
to calculate the entanglement of the superposition, shown in the above equations. 

Fortunately, we have already evaluated the overlap of the two states with the same 
product state |$('C))! as summarized in Eqs. (|36|1 and (!37|l . We only need to compute the 
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sum of the two overlaps weighted by the coefficients cos 6 and sin 6 and maximize the norm 
of the overlap with respect to the parameter C,- As an illustration, we have computed the 
entanglement of superposition for the Ising model with = 10 spins, as shown in Fig. [101 
The resultant entanglement does depend on the detail of the superposition. 

The above exercise can be generalized to all values of r, h, and A^, including the infinite 
limit. If we are interested only in the entanglement per spin, it is now readily easy to see 
from Eq. fl36|) that, in the thermodynamic limit, overlaps for |^o) and |\l'i/2) are identical, 
up to the prefactors f^^ and f^^'^^ . These prefactors do not contribute to the entanglement 
density (although they do matter for the correction in 1/A^), and the entanglement density 
is therefore the same for both |\E'o) and |^I/i/2)- It further follows that, in the thermody- 
namic limit, the results for the entanglement density are insensitive to the replacement of a 
symmetric ground state by a broken-symmetry one. 



E. Finite-size scaling of entanglement at critical points 

Recall that we have performed a finite-size scaling of the entanglement derivative in 
Sec. IVI B 31 in order to extract the critical exponent for the Ising universality class. What 
about the finite-size scaling of the entanglement density itself near criticality? This question 
was recently raised and addressed by Shi and co-workers {ssj , and they have found a scaling 
of the entanglement per site £n as 

^^-^00 + ^ + ^. (66) 

With the formulas of entanglement given above, we can address this question straight- 
forwardly. But as we have discussed earlier, the correction to the entanglement may depend 
on which superposition of degenerate states we use. To avoid the ambiguity, we present 
the finite-size scaling at h = 1 for both states |^E'i/2) (which is the ground state at /i = 1) 
and |\l/o) at various values of r 7^ 0. We shall see that although their £^oo's are practically 
the same, the scaling coefficients b and c are different. One interesting observation is that 
(for r 7^ 0) we obtain b ~ 1.0 for |\I'i/2) and b ^ 0.5 for |\l/o)- The correction b/N has 
recently been shown to related to the Affieck-Ludwig boundary entropy by Stephan and 
54l ]. Furthermore, c is negative for |^I/i/2) and positive for l^&o)- (For r = 



co-workers 



XX model ai h = 1, the entanglement density is identically zero. There is no need to per- 
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0.8 


0.0264164 


1.00001 


-1.69019 


0.0264163 


0.500101 


3.35452 


0.9 


0.0288377 


1.00001 


-1.55332 


0.0288376 


0.50008 


3.08608 


1.0 


0.0311291 


1.00001 


-1.44349 


0.031129 


0.500066 


2.87013 



TABLE I: Table for the scaling parameters b, and c at various r but at fixed h = l for |^'i/2) 
and l^'o)- The results are obtained from fitting the number of spins from 100 to 1000. 

form the finite-size scaling.) To get the scaling fit, we use the size from 100 to 1000, 
and calculate the respective entanglement for fitting Eq. (l66l) . The results are tabulated in 
Table [H Curiously, the c's for 1^1/1/2) are approximately half and negative of those for I^E'o)) 
i.e., c(i/2) ~ -c(o)/2. 

VII. CONCLUDING REMARKS 

In summary, we have quantified the global entanglement of the quantum XY spin chain. 
This model exhibits a rich phase structure, the qualitative features of which are reflected 
by this entanglement measure. Perhaps the most interesting aspect is the divergence in the 
field-derivative of the entanglement as the critical line {h = 1) is crossed. The behavior 
of the divergence is dictated by the universality class of the model. Furthermore, in the 
thermodynamic limit, the entanglement density vanishes on the disorder line (r^ + /i^ = 1). 
The structure of the entanglement surface, as a function of the parameters of the model (the 
magnetic field h and the coupling anisotropy r), is surprisingly rich. We have also discussed 
issues of superposition of entanglement, its independence on the spontaneous breaking of 
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the Z2 symmetry, and finite-size scaling at critical points 

How does ground-state entanglement reflect quantum phase transitions? Can we have 
a generic understanding of the connection? Near quantum critical points, the correlation 
length generally scales as 

i{K)^\h-h\{-\ (67) 

with V denoting the so-called correlation- length exponent. For any density functions, such 
as free energy density and the entanglement density, we expect scaling for the singular part 



behaves as (e.g. Ref. |55l |) 



m ~ (e/a)"'(ci + C2 log(e/a) + ...), (68) 

where d is the dimension of the system, a is the smallest length scale, c's are some constants 
which may depend on which side of critical point we are at, and we have included a possible 
logarithmic correction in its simplest form. The idea is that near criticality the only length 
scale is the correlation length, and any density quantity should scale as 1/volume, i.e., 
l/(^/a)'^, with a being a unit of length. We then expect that in general 



dh 



~ ±a'^\h - hc\'^''-\di + C2 log \h - h,\). (69) 



For the one-dimensional {d = 1) Ising universality class, v = 1, and there is a logarithmic 
divergence, which is consistent with Eq. (1601) . For the XX universality near /i = 1, there is 
no logarithmic correction (i.e., C2 = C2 = 0), the correlation length exponent z/ = 1/2. The 
resultant divergence from the above scaling form is consistent with what we have obtained 
in Eq. (15^ . noting that the coefficient Ci = for /i > /i^. 

We close by pointing towards a deeper connection between the global measure of entan- 
glement and the correlations among quantum fluctuations. The maximal overlap ([1]) can be 
decomposed in terms of correlation functions: 

^max = ^ + ^ H=i {^^' + \ X^^^' ^1 ® '^j) + ■ ■ ■ }' 

i=2 

where translational invariance is assumed and the Cartesian coordinates of r can be taken to 
be (sin^, 0, cos^). The two-point correlations appearing in the decomposition are related to 
a bi-partite measure of entanglement, namely, the concurrence, which shows similar singular 
behavior J to Eq. dSUD- 
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In fact, it was elaborated in Ref. 35| that a singularity of the entanglement can come 
from two types of sources: (i) correlation functions, C^^j'jJ"'"^ = {'^'^]'^[j]^[k] • • • ) for the ground 
state and (ii) parameters r *W, which denote the vectors that maximize the overlap. For 
transverse Ising, which has a standard second-order quantum critical point, (i) correlation 
functions C's are singular but (ii) optimal parameters r*'s are not singular and this explains 
the similar behavior between the GE and the concurrence. It can happen that the singularity 



of the entanglement near a transition can come from different types. In Ref. [35|, it was 
realized that near the so-called Kosterliz-Thouless point (the isotropic antiferromagnetic 
point) of the XXZ model, the correlation length is finite, but (ii) r*'s are singular. It is this 
second point the one that detects non-analyticity in the wavefunction across the transition. 
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Appendix A: Finite-size scaling 



The discussion here follows Barber 50|]. In the vicinity of the bulk critical temperature Tc 
the behavior of a system should depend on y = L/C,{T), where ^(T) is the bulk correlation 
length and L is the characteristic length of the system. How does the divergence of certain 
thermodynamic quantities emerge as the system size L grows? 

1. Algebraic divergence 

Assume that some thermodynamic quantity at L — oo diverges as t = (T—Tc)/Tc oo: 

Poo{T)^C^t-P. (Al) 
Finite-size scaling hypothesis asserts that for finite L and T near Tc, 

Pl{T) ~ rQp{l^'''t), l^oo, i^O, (A2) 
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where I = L/a {a is some microscopic length), i = [T — Tc{L)]. The exponent u can be 
determined by the requirement that the Pl{T) reproduces Poo{T) as Z — > oo. Thus, 

Qp{x) ~ C^x"^, X — >■ oo, (A3) 

and ui — p/v. We consider the case that the finite system does not exhibit a true transition, 
then 

Qp{x)^Qo, x^O. (A4) 

From this we have that at the peak or rounding temperature T^{1) (where Pl reaches the 
maximum or deviates significantly from the bulk value) 

Pz.(r* (/)) - go/''^ /^oo. (A5) 

This means that the behavior of a thermodynamic quantity varies with the system size is 
determined by the bulk critical exponent. 

2. Logarithmic divergence 

Now assume the thermodynamic quantity P{T) diverges logarithmically as 

Poo(T) ~Coolnt, t^O, (A6) 

as in the field-derivative of the entanglement density for anisotropic XY spin chains. The 
finite-size scaling hypothesis in this case is to assume 

Pl{T) - Pl{To) ^ Qpil'^-'i) - Qp{l"%), (A7) 

where Tq is some non-critical temperature and to = (^o ~ Tc{L))/Tc- The hypothesis has 
to recover the / oo limit at fixed T, which requires 

Qp{x) Coo^^x, X ^ OO. (A8) 

Thus in the hmit f — > at fixed large we have 

P^{Tc{L))^-^lnl + 0{l), (A9) 
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if Qp{x) — 0{1) as X — > 0. This allows us to obtain the exponent u by analyzing how 
the divergence develops as the system size I (which is N in our spin-chain entanglement 
problem) increases. 
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